Thermodynamics of a Bose Einstein condensate with free magnetization 
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We study thermodynamic properties of a gas of spin 3 Cr atoms across Bose Einstein conden- 
sation. Magnetization is free, due to dipole-dipole interactions (DDIs). We show that the critical 
temperature for condensation is lowered at extremely low magnetic fields, when the spin degree of 
freedom is thermally activated. The depolarized gas condenses in only one spin component, unless 
the magnetic field is set below a critical value, below which a non ferromagnetic phase is favored. 
Finally we present a spin thermometry efScient even below the degeneracy temperature. 
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A multi-component Bose Einstein condensate (BEG) 
(i.e. a spinor BEG) reveals ground and excited-states 
properties very different from scalar BEGs [H, 0] • Opti- 
cal traps allow the study of such systems as all Zeeman 
states are trapped with almost the same trapping po- 
tential. New quantum phases arise due to an interplay 
between spin dependent contact interactions, linear and 
quadratic Zeeman effects 0]. While these phases have 
been studied in a pioneering set of experiments the 
properties at non zero temperature remain mostly un- 
explored (see however In addition, DDIs, up to 
now neglected in this context may lead to new quantum 
phases and non-trivial spin textures Q. 

In this paper, we study the thermodynamic properties 
of a multi-component gas, made of spin S—3 chromium 
(Gr) atoms, at low temperatures. Our system strongly 
differs from alkali gases due to strong DDIs. One impor- 
tant consequence in the context of spinor gases is that 
magnetization is free 0|- We operate at extremely low 
magnetic fields {qjUbB « fc^T, with /i^ the Bohr mag- 
neton, ks the Boltzmann constant, T the temperature 
and gj = 2 for Gr) which results in non-zero population 
in excited Zeeman states at thermal equilibrium. We 
therefore investigate the phase diagram of spinor bosons 
with free magnetization. 

Our main results are summarized in Figure [T] At large 
temperatures, we observe a normal (A) phase. For large 
enough magnetic fields, and below a magnetization de- 
pendant critical temperature Td , a one component BEG 
(B phase) is obtained. Below a critical magnetic field 
Be [8], corresponding to the field below which spin de- 
pendent contact interactions favor a non ferromagnetic 
ground state [l^l , we observe a multi-component BEG 
(G phase). Figure [1] maps these results onto the phase 
diagram of non interacting spinor bosons as a function of 
temperature and magnetization [ll| . 

The existence of the spin degree of freedom modifies 
statistics and therefore the thermodynamic properties of 
the (25* -I- 1) states come into play. Theoreti- 



1.4- 



1.2- 



1.0- 



t-H 0.8- 
0.6- 



0.4- 



0.2- 



0.0- 




0.0 -0.5 -1.0 -1.5 -2.0 -2.5 -3.0 

Magnetization 

FIG. 1: (Color online) Phase diagram of a spin 3 BEC. The solid lines 
delimitate the 3 phases predicted for a non interacting gas of bosons. A 
phase: thermal gas in each Zeeman component; B phase: BEC only in 
ms — —3; C phase: BEC in all Zeeman components. The histograms 
represent typical experimental population distributions. Black trian- 
gles arc experimental Td measurements (see text). The dashed lines 
apply to regimes below Be (left) or above Be (right). 



feet [13| . They assume that collisions are fast enough for 
thermal equilibrium to be reached, but that interactions 
are weak enough that they do not modify the phase dia- 
gram. These studies generalize results obtained for scalar 
bosons, by introducing a Bose occupation factor for the 
non-condensed particles in each Zeeman state, leading to 
the following dependence for the thermal populations: 
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cal studies of non interacting spinor bosons [ll|, |12| have 
been performed in presence of a purely linear Zeeman ef- 



with LUi the trap frequencies, 1//3 = ksT, and fi the chem- 
ical potential. B is either the experimental magnetic field 
if the magnetization is considered free, or, alternatively, 
an effective magnetic field self-consistently adapted to fix 
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the magnetization thus behaving as a magnetic chemical 
potential. 

Compared to the scalar situation, the general trend 
is a reduction of the critical temperature for condensa- 
tion due to the increased number of degrees of freedom. 
The critical temperature of a polarized sample, Tco, is 
reduced by a factor {2S + ly/^ ( in a 3 dimensional har- 
monic trap) when the magnetization M reaches 0. Two 
different regimes have been studied. When M is fixed 
a double phase transition occurs: below a first criti- 
cal temperature Tci{M) a polarized BEC emerges (phase 
B) i.e. only the lowest energy Zeeman component under- 
goes condensation, while all other components condense 
below a second critical temperature Tc2{M) (phase C). 
This double phase transition is represented by a vertical 
line in the phase diagram of Fig[TJ For a free magnetiza- 
tion [13] , BEC of a non-interacting gas only occurs in the 
lowest energy single particle state {ms — —3 for ^^Cr), 
which corresponds to the B phase: the BEC is always fer- 
romagnetic, at any finite magnetic field, and the C phase 
is always avoided. In this regime the system follows the 
dashed line of Fig [U the function M{T) being calculated 
for a given magnetic field (here 1 mG), which maps the 
results of onto the phase diagram of [lH . 

There are open questions regarding how such theoreti- 
cal frameworks can describe ^^Cr BECs at low magnetic 
fields: the magnetization is free, but the BEC is not 
ferromagnetic below i?c, as we have observed recently 
[lO| . In this paper we show that the theory without in- 
teractions correctly describes thermodynamic properties 
of ^^Cr gases, provided magnetization is considered free 
above i?c, and fixed below. 

To perform our experimental study, we use forced evap- 
oration to prepare a sample of typically 10 to 30 thou- 
sand ^^Cr atoms, at temperatures ranging from 1.5 /^K to 
50 nK. These samples are produced in a crossed optical 
dipole trap [lj|. We repeat the procedure for different 



values of the final energy cutoff set by the evaporative 
ramp, thus varying T; then the trap is recompressed at 
final trap frequencies independent of T . We control the 
B field at the 100 /xG level, and the magnetization of 
the cloud, M = J2'i=-3Pi ^ * (where pi is the relative 
population of the rrii Zeeman state) is measured by a 
Stern-Gerlach procedure [l3|- By varying the value of B 
and T, we can investigate, after thermal equilibrium is 
reached, properties of a spin 3 gas for different magneti- 
zations and temperatures. 

For 6 different values oi B, we have measured the 
steady state magnetization M as a function of T, and 
the condensed fraction. We deduce six values of Td (for 
which the condensed fraction departs from zero), which 
are reported in Fig[T] The observed decrease of Td as M 
increases from its minimal value —3 corresponds to a re- 
lease of the spin degrees of freedom: the number of atoms 
in the most populated ms = —3 state is lower than the 
total number of atoms, so that saturation is achieved at 
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FIG. 2: (Color online) Evolution of the magnetization with tempera- 
ture. Triangles: above Be {B — 0.9 mG), the magnetization M slowly 
departs from zero as the temperature is lowered. When condensation is 
reached M decreases faster, and reaches about —3 for the lowest tem- 
peratures; the data agree with a theory without interactions and free 
magnetization (solid line). Diamonds: below Be (here B — Bmin): the 
magnetization remains almost constant and significantly different from 
—3, even at extremely low temperature (the horizontal solid line shows 
its average value). 



lower temperature. Fig [T] shows good agreement between 
our results for Tci{M) and the theory of non interacting 
spinor bosons. Interestingly, the agreement between this 
simple model and the experimental results is good, while 
we explore values of B both above and below Be, corre- 
sponding to BEC of very different nature at T = 0. This 
shows that Td is relatively insensitive to interactions for 
our system, in both B regimes; this is expected in the 
case of 5 = 1 atoms [11]. 

Below Td, we find two entirely different behaviors of 
the system depending on the relative value of B and Sc- 
This is illustrated in Fig [21 where we show experimental 
values of the gas magnetization as a function of tem- 
perature for two different magnetic fields. Above B^, 
when temperature is decreased, the magnetization of the 
purely thermal (for T > 400 nK) gas first progressively 
decreases from about zero (all the spin components being 
almost equally populated), to negative values, as the neg- 
ative spin components become more populated. Below a 
certain temperature (w 400 nK), we observe a kink in 
magnetization as well as a bimodal spin population dis- 
tribution (see the experimental histogram characteristic 
of the B phase in Fig[T]). At this temperature, a BEC is 
reached, as independently shown by the appearance of a 
bimodal momentum distribution measured after a time of 
flight. We interpret the kink in magnetization by the fact 
that only atoms in the ms = —3 component condense, 
so that magnetization decreases faster with temperature. 
M = —3 is reached for the lowest temperatures, when 
almost all the atoms are Bose condensed. Fig [5] shows 
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the results of the theory without interactions and with 
free magnetization, and the good agreement with our ex- 
perimental data. We hence pinpoint a spontaneous po- 
larization of the cloud as condensation is reached, due 
to its ferromagnetic nature. Our data only reveal the 
equilibrium properties and it would be most interesting 
to study the spontaneous polarization of the condensed 
part due to DDIs if the temperature is lowered below Td 
fast compared to magnetization dynamics. 

Below Be on the other hand, the results strongly de- 
part from predictions of thermodynamics of non inter- 
acting bosons with free magnetization. We show in Fig 
[2] results for the lowest magnetic field we can achieve 
{B = Bmin < 0.1 mG). The magnetization does not 
reach —3, even for the lowest temperatures, whereas the 
condensate fraction reaches almost 1 in this situation (see 
Fig|3]). Magnetization remains in fact almost constant, 
allowing to enter a C-like phase at low temperatures. The 
population distribution at the lowest temperature is rep- 
resented by the lower histogram in FigHJ For these data, 
we have checked that the thermal populations are neg- 
ligible. It is worth noting that even though we enter 
a spinor phase, the spin distribution that we obtain is 
not the one predicted for the absolute ground state of 
chromium Q. This may indicate that thermal equilib- 
rium is not reached in this specific case, and that the 
BEG is in a metastable spin state. 

A special emphasis is made for the temperature depen- 
dance of the condensed fraction rj below Be in FigO The 
filled triangles correspond to data taken at large magnetic 
field {B = 20 mG), while the open triangles correspond 
to B = Bjnin < Be- At large B field, the cloud remains 
polarized. In this case we measure a critical temperature 
Tco, close to the value predicted for a non-interacting gas, 
as expected for weak interactions [lil - flTj . For B < Be 
we observe a dramatic change: t] is significantly smaller, 
at any temperature below the critical temperature Tcq, 
than at high field. Large values of rj are only recovered 
for the lowest temperatures. 

To explain these low B fields results, we compare them 
with a theory with no interactions, but considering a 
fixed magnetization, as is empirically almost the case be- 
low Be (see Fig [2]). Gorresponding results are shown 
by the solid line in FigJSl illustrating the double phase 
transition discussed in [ill ]. In this experimental config- 
uration, we measure Td ~ 300 uK. Following the theory 
in ll|, all Zeeman components should condense below 
Tc2 « 150 nK, and although we do not measure Tc2 we do 
find a multi-component BEG at the lowest temperatures 
[lo| . The agreement between experimental data and the 
predictions of the model in Fig [3] is satisfactory, hence 
showing a hint for the double phase transition. Thermo- 
dynamics of non interacting bosons with (almost) fixed 
magnetization hence seems to correctly describe our sys- 
tem. We however emphasize that a spinor G phase is 
only achievable below Be because then the ground state 




0.2 0.3 0.4 
Temperature (|iK) 

FIG. 3: (Color online) Evolution of the condensed fraction r] with 
temperature. Filled triangles: for a polarized sample {B — 20 mG), 
the usual 1 — (T/Tco)"^ (corresponding to the right solid line) is ob- 
tained. Open triangles: at Bmin (below Be), the condensed fraction 
gets smaller at every temperature, as the spin degree of freedom is un- 
frozen. The left solid line corresponds to a theory with no interactions 
but fixed magnetization. 



is not ferromagnetic, due to spin dependent contact in- 
teractions [10]. 

Finally, we take advantage of results obtained for B 
fields just above Be (i.e. in the B phase) in order to 
investigate a spin thermometry, efficient down to tem- 
peratures much smaller than Teo- As explained above, 
for B > Be the BEG is only in ms = —3, while the other 
7715 components, purely thermal, obey a Bose statistics. 
A small value of B allows thermal population in all 7775 
components. We show in the inset of Fig |4] the different 
populations measured for B = 0.5 mG, at a temperature 
of 160 nK. The Bose statistics favors the negative val- 
ues of 7775, while the large population in —3 emphasizes 
the ferromagnetic nature of the BEG. We derive the spin 
temperature, T^pm, by fitting the ms populations (for 
7775 > —3) to an exponential [l8| . 

In Fig |4] we compare Tspm to the temperature Ttof 
measured by a standard bimodal fit after a time of flight 
of 5 ms (without Stern-Gerlach separation). We observe 
a good general agreement between the two temperatures, 
indicating that spin and spatial degrees of freedom are at 
equilibrium. The error bars in Tspin are a combination of 
the effect of the uncertainties on B (dominant at low T) , 
and of the statistical errors on the flt of the thermal pop- 
ulations [1^. At T way below Teo, it is difficult to extract 
accurate ternperatures from a bimodal fit of the velocity 
distribution [20| : the thermal cloud barely separates from 
the condensed part for /cbT < /i/3, hence the increase in 
the error bars of Ttof at the lowest temperatures. At 
low T, Tspin provides a more accurate measurement than 
Ttof, showing the interest of this method. 

This spin thermometry, relevant for dipolar gases, com- 
pares well with other techniques, such as noise thermom- 

[23. It bears 



etry with BEGs coupled in a double well 
similarities with spin thermometry using a magnetic field 
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FIG. 4: Comparison of the cloud temperatures deduced from a time 
of flight analysis (Ttof) with those of a spin components analysis 
(Tspin), for B — 0.9 mG. Inset: spin components populations. The 
magnetic field is set above Be (B — 0.5 mG), so that only the thermal 
fraction gets depolarized. The temperature is below Td (T — 160 
nK), hence the large ms — —3 component. The other ms components 
almost obey a Boltzmann statistics I18II . as shown by the LogPlot, from 
which Tg^in is deduced. 



gradient [22|, but avoids the use of a gradient. As the 
minimal measurable temperature scales with i?c, this 
thermometry will therefore be most interesting for ei- 
ther ferromagnetic atoms, or atoms with small Be] nev- 
ertheless, interactions between thermal atoms and the 
BEC could introduce systematics effects. Besides, one 
may take advantage of the selective depolarization of the 
thermal gas in the B phase to develop a cooling method, 
based on a configuration where only the mg = — 3 com- 
ponent is trapped: this would provide a selective loss pro- 
cess for the thermal cloud, which should result in cooling. 

In conclusion, by working at low temperature and low 
magnetic field we have measured unique thermodynamic 
features of a gas of atoms with non zero spin, and free 
magnetization. We have evidenced strong differences in 
the behavior of the gas above or below Be, the critical 
field separating ferromagnetic and unpolarized phases. 
We have mapped our results onto the phase diagram of 
a non interacting S=3 Bose gas. Below Be a metastable 
spinor C-like phase is observed. Above Be, the BEC is 
always polarized, which we use to introduce a thermom- 
etry accurate even below the critical temperature. 
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